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摘 要 : 选取 2000—2018 “4 MODIS-NDVI 数据 ,采用 Sen + Mann-Kendall # # 4 #7. 变异 系数 法 、 
Hurst 指数 、 偏 相关 分 析 及 残 差 分 析 等 方法 ,分 析 塔 里 木 河 流域 植被 覆盖 时 空格 局 及 气候 因子 和 人 


类 活动 对 植被 变化 的 影响 。 结 果 表 明 :(1) HE 


木 河 流域 植被 覆盖 总 体 呈 明显 增加 趋势 ,以 2008 


年 为 转折 存在 明显 阶段 变化 。(2) 植被 分 布 存在 明显 地 域 差异 ， 北 高 南 低 , 西 高 东 低 ” ,植被 覆盖 


kt 


区 分 布 于 山体 多 的 地 带 、 绿 洲 及 绿洲 荒漠 交错 带 。(3) 全 区 植被 覆盖 趋势 以 基本 不 变 居多 。 山 


体 等 植被 覆盖 高 区 ,植被 活动 响应 显著 。 全 区 大 部 分 区 域 将 保持 现 有 的 稳定 趋势 。(4) 气温 和 降 


水 量 对 NDVI 的 综合 影响 / 


东北 向 西南 逐渐 增强 , 且 降 水 量 对 NDVI 的 影响 更 为 明显 。 人 类 活动 对 


植被 变化 产生 积极 影响 的 区 域 主要 分 布 于 绿洲 、 绿 洲 蕊 漠 过 渡 带 以 及 塔里木 河 下 游 附 近 。 
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植被 作为 陆地 生态 系统 重要 的 主体 之 一 ,与 大 


气 圈 .土壤 圈 、 水 圈 的 物质 循环 和 能 量 流动 紧密 关 
联 '"1。 植 被 覆盖 变化 改变 地 表土 壤 下 垫 面 属性 ,对 
水 土 保持 .气候 调节 以 及 生态 环境 的 稳定 都 起 到 关 
键 作用 ,植被 覆盖 度 同 植被 繁 成 程度 和 植被 进行 光 
合作 用 的 面积 有 明显 的 对 应 关系 '”。 归 一 化 植被 
指数 (Normal Difference Vegetation Index , NDVI) #2: H 
前 常用 的 表征 植被 覆盖 变化 的 指标 ,能够 揭示 区 域 
环境 状况 的 演化 与 变迁 ”。 将 NDVI 数据 应 用 于 监 
测 植被 变化 特征 ,了 解 植被 时 空 覆盖 格局 .过程 , 进 
而 为 探索 提高 区 域 生 态 环境 质量 的 方法 .路 径 ,改善 
人 们 居住 环境 .促进 地 域 社 会 经 济 发 展 提供 参 考 。 
长 期 以 来 由 于 气候 变化 和 不 合理 的 开发 利用 ， 
塔里木 河 下 游 河 道 断 流 ,绿色 走廊 不 断 误 退 EVE AG 
加 剧 ,严峻 的 生态 环境 问题 严重 制约 了 地 域 社会 经 
济 发 展 ”。2001 年 起 ,塔里木 河流 域 管理 局 先后 实 
施 输 水 工程 灌区 节 水 .河道 治理 等 一 系列 生态 综合 
治理 工程 , 旨 在 缓解 塔里木 河流 域 生 态 环境 日 趋 恶 
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化 。 利 用 遥感 数据 对 塔里木 河流 域 植被 覆盖 状况 的 
研究 已 有 许多 成 果 : 郭 继 凯 等 “基于 2001—2013 年 
MODIS-NDVI 数据 和 塔里木 河流 域 气象 数据 ,利用 
一 元 线性 回归 法 、 残 差分 析 法 等 方法 ,分 析 了 塔里木 
河流 域 植被 覆盖 变化 趋势 .气候 变化 以 及 人 类 活动 
对 植被 改善 区 和 植被 退化 区 的 影响 ; 罗 敏 等 ' 基于 
2000—2013 年 MODIS-NDVI 数据 和 塔里木 河流 域 
同期 气候 数据 ,使 用 一 元 线性 回归 分 析 、 偏 相关 分 析 
以 及 残 差 分 析 等 方法 ,研究 了 塔里木 河流 域 生长 季 
NDVI 的 时 空 变化 特征 及 其 影响 因素 分 析 ; 管 文 町 
等 中 基 于 2007—2017 年 MODIS-EVI 数据 ,运用 线 
性 回归 和 方差 分 析 等 方法 ,对 塔里木 河流 域 植被 的 
变化 规律 进行 动态 监测 及 趋势 分 析 ; 陈海燕 等 咏 基 
于 1977—2014 年 近 40 a 的 5 个 不 同时 期 的 Landsat 
数据 ,采用 NDVI 计算 模型 一 元 线性 回归 分 析 等 方 
法 ,分 析 了 塔里木 河 干流 荒漠 河岸 林 植 被 带 的 年 际 
变化 及 空间 格局 态势 ;高 庆 等 ”基于 1986 一 2017 年 
Landsat 数据 ,利用 波动 分 析 ,趋势 分 析 等 方法 ,研究 
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了 塔里木 河 下 游 区 域 植被 时 空 变化 特征 ,并 对 其 影 
响 因素 进行 分 析 ; 闫 俊杰 等 '” 基于 2000—2014 年 
MODIS-NDVI 数据 ,使 用 Sen + Mann-Kendall 趋 热 分 
析 等 方法 ,对 塔里木 河流 域 植被 年 际 变化 特征 及 空 
间 趋 势 进行 分 析 。 

上 述 研究 成 果 主 要 在 不 考虑 其 他 线性 因素 影响 
下 ,研究 气温 、 降 水 量 和 人 类 活动 对 植被 NDVI 的 影 
响 ,实际 上 ,植被 的 生长 变化 受 诸 多 因素 综合 影响 ， 
单一 的 某 种 因素 往往 不 能 解释 植被 复杂 的 变化 过 
程 。 再 者 ,很 少 研究 者 对 塔里木 河流 域 植 被 覆盖 空 
间 格 局 态势 进行 综合 分 析 。 鉴 于 上 述 情况 ,本文 基 
于 MODIS-NDVI 数据 ,利用 偏向 关 分 析 方 法 研究 气 
温和 降水 量 对 塔里木 河流 域 植被 变化 的 综合 影响 ， 
利用 残 差分 析 方 法 分 析 人 类 活动 对 植被 变化 的 影 
响 ,进一步 反映 NDVI 残 差 在 人 类 活动 相对 较 弱 区 域 
还 受 其 他 线性 变化 因素 影响 ,结合 Sen + Mann-Ken- 
dal 趋势 分 析 、 变 异 系数 法 Hurst 指数 等 方法 分 析 塔 
里 木 河 流域 植被 覆盖 时 空格 局 变化 特征 及 趋势 。 


1 数据 与 方法 


1.1 研究 区 概况 

塔里木 河流 域 位 于 新 疆 南部 , 环 塔里木 盆地 与 
山地 相 接 , 北 依 天 山 , 西 临 帕 米尔 高 原 , 南 靠 昆仑 山 、 
阿尔 金山 , 东 至 库 鲁 克 塔 格 山 ,塔克拉玛干 沙漠 位 于 
流域 中 部 '" 。 受 水 资源 开发 利用 影响 ,目前 逐渐 形 
成 流域 的 阿克苏 河 、 和 田 河 、 叶 尔 羌 河 \ 开 都 河 4 条 
源流 和 塔里木 河 干 流 的 格局 ,塔里木 盆地 边缘 与 流 
域 之 间 形 成 冲积 扇 和 冲击 平原 ,并 有 大 片 绿洲 分 
WU 。 塔 里 木 河 流域 深 居 内 陆 , 远离 海洋 ,高山 环 
绕 , 属 于 温带 干旱 大 陆 性 气候 , 光 热 资 源 丰 富 ,流域 
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内 干旱 少雨 ,多 年 年 均 降水 量 不 超过 120 mm ,蒸发 
强烈 , 且 年 .日 温差 大 ,年 均 气 温 除 高 寒山 区 外 多 在 
3.2 ~11.5 C, WEARER! 。 

塔里木 河流 域 特 殊 的 地 理 区 位 使 得 水 资源 匮 
Z ,植被 极度 稀 琉 且 覆 盖 低 ,戈壁 .沙漠 裸 地 面积 大 ， 
具备 典型 的 干旱 区 生态 环境 特征 王 。 塔 里 木 河流 
域 整体 按 行政 区 可 划分 为 7 大 植被 覆盖 区 :哈密 地 
区 吐鲁番 地 区 ` 巴 音 郭 楞 蒙 古 自 治 州 ( 简 称 巴 州 )、 
阿 元 苏 地 区 、 克 孜 勤 苏 柯尔克孜 自治 州 ( 简 称 克 
州 ) 喀什 地 区 .和 田地 区 (图 1)。 
1.2 研究 数据 

本 研究 中 遥感 数据 由 美国 地 球 观测 系统 数据 和 
信息 系统 EOSDIS (Earth Observing System Data and 
Information System) 提供 的 2000—2018 年 MODIS/ 
MOD13Q1 植被 指数 产品 ,16 d 最 大 值 合成 .250 m 
空间 分 辨 率 。 利 用 MRT 软件 从 MOD13Q1 产品 中 
提取 NDVI 数据 ,经 灸 内 .裁剪 MVC 合成 等 预 处 理 ， 
得 到 研究 区 NDVI 数据 集 。 以 % 月 作为 研究 时 段 ， 
此 时 段 区 域 NDV 达到 一 年 中 最 大 , 且 影像 云 污染 
较 少 ,能 够 连续 对 研究 区 植被 时 空 变 化 特征 进行 分 
析 。 气 候 数 据 来 源 于 中 国 气象 数据 共享 中 心 , 包 括 
研究 区 19 个 气象 站 点 2000 一 2018 年 同期 气温 及 降 
水 数据 。 
1.3 研究 方法 
1.3.1 Sen + Mann-Kendall 趋势 分 析 
Mann-Kendall 方法 是 Sen 趋 热度 与 Mann-Kendall 检 
验 结合 ,以 增强 长 时 序 趋势 分 析 抗 噪 性 ,提高 检验 结 
REME o Sen 趋势 度 计 算 公式 为 : 


Sen 十 


[NDVI - NDVI, _ 
B = median a Vi>i (1) 
90°E 95°E 100°E 


图 例 
。 站 点 
一 行政 区 划 


一 主要 河流 
NDVI 


wr 高 :0.990 
低 :0.001 


90°E 95°E 100°E 


图 1 研究 区 示意 图 


Fig.1 Location of the study area 


孙 天 瑶 等 :2000 一 2018 年 塔里木 河流 域 植 被 覆盖 时 空格 局 


式 中 :NDVI; 和 NDVI, 为 样本 时 序数 据 集合 ;median 
为 中 值 函 数 ;1 <i <j <n ij 分 别 为 集合 不 同时 间 长 
度 ,n 为 时 序 长 度 ; V 为 任意 符号 。B 用 以 量化 单调 
趋势 , 当 B >0 时 ,反映 这 一 时 间 序 列 呈 单调 上 升 趋 
热 , 反 之 呈 单 调 下 降 趋势 。 

使 用 Mann-Kendall 进行 NDVI 序列 趋势 检验 ， 
设 检验 统计 量 Z, 其 计算 公式 为 : 

S-1 


S>0 
Vvar(S) 
Za @ ga (2) 
_95+l gc0 
~ var( S) 
s S $ sgn( NDVI, ~ NDVI, ) (3) 
‘tere n(n -1)(2n +5) (4) 


18 
1 NDVI, - NDVI, >0 


sgn (NDVI, ~ NDVI,) = NDVI, - NDVI, =0 (5) 


-1 NDVI, - NDVI, <0 


IUP NDVI, 和 NDVI 分 别 为 样本 时 序数 据 集合 ;var 
为 计算 方差 函数 ;sgn 为 符号 函数 ;1 <i<k<n,i,k 
es ae 根据 时 序 长 度 n 值 大 小 
的 不 同 , 当 显 著 性 检验 统计 量 $ 近似 服从 标准 正 态 
分 布 时 ,使 用 Z er 在 给 定 显著 性 
水 平 a 下 ,在 正太 分 布 表 中 查 临 界 值 Zu IZI 
<Z up 时 , 即 趋势 不 显著 ; 当 121 >Z,_.., WAA 
趋势 显著 。 
1.3.2 Hurst 4% Hurst 指数 用 于 定量 描述 长 时 
序数 据 的 可 持续 性 ,判断 时 间 序 列 遵从 随机 游 走 还 
是 存在 持续 性 。 其 计算 原理 如 下 : 
时 间 序 列 | NDVI(t) } ,t=1,2,… 
正 整 数 7 宇 1, 定 义 该 序列 均值 序列 : 


,入 ,对 于 任意 


<NDVI>, = +5 NDVI(1) page CS 
Tif 


累积 离 差 : 
X(t,7) = $ (NDVI(u) - <NDVI>,) ,1<t<7 (7) 
极 差 : 
R(7) = max X(t,7) - min X(1,7) 
7=1,2,-°,N (8) 
标准 差 : 


S(r) = [= vovo - <NDVI>,)*]", 


7T=1,2,，…,N (9) 

若 存 在 R(T)/S(7) x (cr)”,c HARE, H 即 
为 Hurst 指数 。 当 0.00 <H <0. 50 ,表示 时 间 序 列 具 
有 有 反 持 续 性 , 即 未 来 变量 与 过 去 变化 趋势 相反 ,万 越 
小 , 反 持续 性 越 强 ; 当 H =0. 50 时 ,表示 时 间 序 列 为 
相互 独立 的 随机 序列 ; 当 万 >0. 50 ,表示 未 来 变量 后 
过 去 变化 趋势 一 致 ,有 越 大 ,持续 性 越 强 。 
1.3.3 变异 系数 法 ”变异 系数 能 够 衡量 时 间 序 列 
偏离 常态 的 程度 ,用 来 反映 植被 覆盖 波动 变化 。 其 
值 越 大 ,表明 数据 波动 性 越 大 , 越 不 稳定 ;反之 数据 
分 布 较为 集中 ,稳定 性 强 '“ 。 其 计算 公式 为 : 


O NDVI 


NDVI 
BCH: CV yyy) HEAR! NDVI 的 变异 系数 ;0 ww 表示 
时 序 NDVI 标准 差 ;ND 表示 均值 。 
1.3.4 偏 相 关 分 析 ” 偏 相关 分 析 是 指 在 三 个 或 多 
个 变量 的 相互 作用 中 ,控制 其 他 变量 影响 ,研究 其 中 
两 个 变量 的 相关 性 程度 的 过 程 "" 1 。 其 计算 公式 为 : 


CV ypv = 


(10) 


T, 一 六 


ase ao (11) 
Ue JO- 
式 中 :rw Pee nm- 分 别 为 变量 x 和 y、 变 量 x 与 z、 变 量 y 


与 z 间 的 相关 系数 
1.3.5 残 郑 分 析 法 ” 残 差 分 析 利 用 残 差 提供 的 信 
息 ,来 考察 模型 假设 的 合理 性 及 数据 可 靠 性 。 根 据 
气候 要 素 与 NDVI 的 相关 关系 分 离 出 气候 要 素 对 
NDVI 变化 的 影响 , 残 差 变化 趋势 用 来 表示 人 类 活动 
(放牧 农业 ,城市 化 等 ) 和 其 他 线性 变化 因素 的 影 
mg"), BEARN: 

e=y-y' (12) 
式 中 :7y 为 NDVI 实测 值 ;y 为 NDVI 模型 预测 值 ;e 
为 残 差 。 


2 结果 与 分 析 


2.1 塔里木 河流 域 植被 年 际 变化 特征 
2000 一 2018 年 塔里木 河流 域 植 被 覆盖 存在 一 
定 的 波动 ,总 体 呈 现 较 为 明显 的 增加 趋势 ,NDVI 值 
介 于 0. 105 ~ 0. 224 ,年 均 增 长 率 为 4.78% (图 2)。 
植被 覆盖 存在 明显 的 阶段 变化 ,2008 年 以 前 植被 变 
化 由 增加 趋势 转向 下 降 ,之 后 迅速 增长 ,但 在 2014 
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图 2 塔里木 河流 域 NDVI 年 际 变化 曲线 
Fig.2 Annual variation curve of NDVI 


in the Tarim River Basin 


年 有 一 个 较 大 的 下 降 , 不 及 2008 4E NDVI 值 最 低 。 
其 中 克 州 植被 增加 最 为 显著 ,年 均 增 长 率 为 9.42% ， 
其 次 是 哈密 地 区 7. 12% , 巴 州 地 区 3.95% 。 和 田地 
区 植被 增长 率 缓慢 , 这 与 其 毗邻 昆仑 山北 坡 属 暖 温 
带 塔里木 芒 漠 .植被 基数 较 小 有 关 ,但 年 均 增 长 率 亦 
以 2.01% 上 升 ( 表 1)。 


R1 塔里木 河流 域 NDVI 年 际 变化 特征 统计 
Tab.1 Statistical characteristics of annual variability 
of NDVI in the Tarim River Basin 


xi NDVI 年 际 变化 年 均 增长 率 
2000 4Æ 2008 4E 2018 年 / % 
哈密 地 区 0.17 0.27 0.40 7.12 
吐鲁番 地 区 0.18 0.25 0.29 3.22 
巴 州 0.20 0.27 0.35 3.95 
阿克苏 地 区 0.22 0.33 0.35 3.11 
克 州 0.19 0.33 0.53 9.42 
喀什 地 区 0.21 0.29 0.37 4.01 
和 田地 区 0.34 0.32 0.47 2.01 


2.2 塔里木 河流 域 空间 格局 变化 

2.2.1 室 间 分 布 基本 特征 ”塔里木 河流 域 植被 覆 
盖 增 长 显著 (图 3) 。 植 被 覆盖 由 流域 外 侧 向 流域 中 
心 一 塔里木 盆地 腹部 逐年 扩大 , 低 植被 覆盖 区 逐渐 
减少 ;植被 增加 主要 集中 于 库尔勒 BSS .阿克苏 、 库 
车 喀什 地 区 以 及 天 山中 部 , 库 鲁 克 塔 格 山 附 近 植 被 
变化 相对 较 小 ;全 区 植被 状况 相 比 2000 年 均 有 所 提 
高 ,如 同心 年 均 增长 率 由 2008 年 2.11% 上 升 为 
4.56% 。 塔 里 木 河流 域 植被 覆盖 受 区 域 地 貌 格局 影 
响 较 为 明显 ,天山 山脉 、 帕 米尔 高 原 以 及 昆仑 山脉 东 
部 等 附近 形成 明显 的 植被 相对 高 区 ,塔里木 盆地 腹 
部 、 东 部 库 鲁 克 塔 格 沙漠 以 及 阿尔 金山 脉 等 附近 植 
被 覆盖 相对 较 低 ; 受 水 热 条 件 限 制 ,流域 内 植被 覆盖 
北部 高 于 南部 ,西部 优 于 东部 和 东北 部 。 
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图 3 2000—2018 年 塔里木 河流 域 
植被 覆盖 空间 分 布 


Fig.3 Spatial distribution of vegetation coverage 


in the Tarim River Basin from 2000 to 2018 


2.2.2 室 间 格局 趋势 及 稳定 性 ”从 空间 分 布 上 
看 , 越 靠近 山体 ,或 者 山体 越 大 的 地 方 , 植 被 覆盖 增 
加 趋势 越 明显 。 以 塔里木 河 干流 、 开 都 河 、 阿 克 苏 
河 、 克 孜 勒 苏 河 与 和 田 河 形成 明显 的 植被 绿化 区 。 
流域 沿岸 较 流域 中 心 植被 改善 显著 , 呈 改 善 趋势 区 
域 主要 分 布 于 哈密 、 库 尔 勒 .看 背 、 阿 克 苏 、 库 车 、 阿 
图 什 、 乌 恰 喀什 地 区 以 及 天 山南 脉 。 植 被 呈 退 化 趋 
势 区 域 主要 分 布 于 天 山中 部 的 巴 音 布鲁克 草原 .由 
米尔 高 原 以 及 昆仑 山东 段 等 地 区 。 阿 克 苏 地 区 北 
部 .昆仑 山中 段 以 及 库 鲁 克 塔 格 沙漠 等 附近 也 有 人 少 
数 区 域 出 现 植被 覆盖 减少 的 趋势 ,同时 植被 覆盖 改 
善 区 域 中 也 存在 零星 的 退化 区 域 。 塔 里 木 盆 地 腹 
部 、 库 鲁 克 塔 格 沙 漠 以 及 阿尔 金山 等 地 区 植被 覆盖 
变化 处 于 基本 不 变 趋势 。 
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4 2000—2018 年 塔里木 河流 域 NDVI 
空间 趋势 及 稳定 程度 
Fig.4 Spatial trends and stability of NDVI in the Tarim 
River Basin from 2000 to 2018 


R2 2000—2018 年 塔里木 河流 域 NDVI Sen + 
Mann-Kendall 趋势 分 析 统 计 
Tab.2 Statistical results of Sen + Mann-Kendall of 
NDVI in the Tarim River Basin from 2000 to 2018 
NDVI 变化 趋势 程度 比重 / % 
B=0.000 5,121>=1.96 明显 改善 2. 13 
B=0.000 5,171 <1. 96 轻微 改善 35. 84 
-0. 000 5 <8 <0. 000 5 ,1Z1 <1. 96 基本 不 变 50. 80 
B<0.000 5,1Z1 <1. 96 轻微 退化 10. 91 
BO. 000 5,127] =1. 96 严重 退化 0.33 


塔里木 河流 域 植被 覆盖 趋势 以 基本 不 变 区 域 居 
多 (50. 80% ) , 植被 覆盖 呈 改 善 趋势 区 域 比重 为 
37.97% ,轻微 和 严重 退化 趋势 区 域 比重 为 11.24% , 
从 总 体 上 看 ,改善 趋势 区 域 主要 集中 于 流域 内 附近 
水 资源 比较 丰富 的 绿洲 ,退化 趋势 区 域 主要 分 布 于 
巴 音 布鲁克 草原 和 昆仑 山东 段 等 植被 基数 相对 较 
大 ,植被 生长 变化 响应 显著 的 地 区 。 

塔里木 河流 域 植被 覆盖 稳定 性 不 同 地 域 存在 明 
显 差异 (图 4b) ,总 体 上 处 于 “流域 外 侧 高 于 内 侧 ,高 
低 波动 并 存 , 中 度 波动 居多 ”的 态势 。 波 动 程度 高 


绿洲 ,和 田 绿 洲 等 地 区 ,与 植被 改善 区 基本 吻合 ; 
等 波动 变化 区 域 多 集中 分 布 于 绿洲 与 荒漠 交错 带 ; 
相对 低 波动 变化 以 下 区 域 主要 分 布 于 塔克拉玛干 沙 
漠 、 库 鲁 克 塔 格 沙漠 等 人 口 稀 少 . 水 资源 稀少 其 
漠 区 。 

塔里木 河流 域 植被 波动 程度 高 区 域 比 重 为 
40.02% ,中 等 波动 变化 以 下 区 域 比重 为 59. 98% , 
综合 可 以 发 现 ,水 分 条 件 是 控制 流域 内 植被 生长 变 
化 的 主要 因素 :山区 以 及 绿洲 高 波动 区 ,多 以 林地 覆 
盖 , 水 分 条 件 充足 ,气温 趋 于 变 暖 ,植被 生长 变化 响 
应 显著 ;绿洲 荒漠 交错 带 , 受 水 热 条 件 及 车 漠 气 候 限 
制 ,植被 波动 相对 不 高 ;人 口 稀少 ,水 资源 稀少 的 其 
漠 区 ,植被 变化 相对 低 , 其 大 幅 波动 变化 可 能 性 
较 小 。 
2.2.3 植被 覆盖 可 持续 性 分 析 研究 区 NDVI 的 
Hurst 指数 均值 为 0. 65 ,表明 塔里木 河流 域 植被 变 
化 将 继续 保持 目前 变化 趋势 。 从 2000—2018 年 , 反 
持续 序列 比重 为 30. 83% , 主要 分 布 于 巴 音 布 鲁 克 
草原 . 帕 米 尔 高 原 . 昆 仑 山东 段 、 塔里木 河上 中 游 附 
近 ,未 来 可 能 出 现 退 化 趋势 ,这 可 能 由 于 山区 和 水 系 
附近 植被 覆盖 相对 高 ,变化 波动 大 ,植被 活动 稳定 性 
低 。 持 续 性 序列 比重 为 69. 17% , 主要 集中 于 流域 
绿洲 地 带 一 库尔勒 绿洲 、 渭 干 河 三 角 洲 .阿克苏 河 冲 
击 平原 绿洲 \ 喀 什 三 角 洲 绿洲 、 叶 尔 羌 河 绿洲 、 和 田 
绿洲 ,绿洲 若 漠 交错 带 以 及 荒漠 区 (图 5)。 

塔里木 河流 域 呈 正 向 持续 序列 几乎 覆盖 了 整个 
区 域 ,加 之 2000—2018 年 植被 呈 退 化 趋势 比重 仅 占 
11.24% ,由 此 可 预测 塔里木 河流 域 大 部 分 区 域 植被 
将 持续 改善 。 弱 持续 性 序列 (0.35 ~0.65) 比重 为 
63.17% ,几乎 占 一 半 , 表 明 植 被 变化 持续 程度 较 弱 。 
这 种 弱 持 续 性 进一步 说 明 流域 内 植被 覆盖 情况 是 由 
自然 条 件 和 人 类 活动 共同 影响 的 ,单纯 自然 或 人 类 
干预 修复 ,其 恢复 序列 应 具有 较 强 持续 性 ,如 自然 维 
持 的 森林 系统 ,人 类 维持 的 农田 系统 。 
2.3 塔里木 河流 域 气候 因子 对 植被 变化 的 影响 

对 塔里木 河流 域 19 个 气象 站 的 气温 与 降水 数 
据 进行 线性 拟 合 ( 图 6)。19 a 来 ,塔里木 河流 域 降 
水 量 呈 明显 上 升 趋势 ,气温 呈 微 弱 上 升 趋 势 ,气候 环 
境 正在 逐步 由 暖 干 向 暖 湿 化 趋势 转变 。2000 一 2008 
年 ,塔里木 河流 域 平 均 气 温 呈 缓慢 上 升 趋势 ,年 降水 


区 域 主要 分 布 于 帕 米 尔 高 原 、 天 山南 脉 \, 昆 仑 山中 东 
段 及 流域 绿洲 区 一 库尔勒 绿洲 , 渭 干 河 三 角 洲 Ba 
苏 河 冲击 平原 绿洲 喀什 三 角 洲 绿洲 \. 叶 尔 羌 河 平 原 


量 呈 明显 下 降 趋势 ;2009 一 2018 年 ,塔里木 河流 气 
温 依旧 呈 上 升 趋势 ,降水 量 呈 明显 上 升 趋势 。 
为 了 进一步 分 析 气 候 因 子 对 植被 变化 的 影响 ， 
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Fig.5 Spatial distribution of Hurst Index in the Tarim River Basin from 2000 to 2018 
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图 6 2000—2018 年 塔里木 河 流域 气温 
和 降水 量变 化 趋势 


Fig.6 Trend of temperature and precipitation in the Tarim 


River Basin from 2000 to 2018 


采用 相关 分 析 法 分 析 气 温 以 及 降水 量 因子 对 塔里木 
河流 域 NDVI 的 影响 。2000 一 2018 年 塔里木 河流 域 
气温 与 WD 人 的 相关 系数 ,与 降水 量 与 NDVI 的 相关 
系数 ,区 域 差异 明显 (图 7)。 塔里木 河流 域 NDVI 
在 大 部 分 区 域 (72.92% ) 与 气温 呈 明 显 的 负 相 关 关 
系 , 主 要 分 布 在 吐鲁番 地 区 阿克苏 地 、 克 州 .喀什 地 
区 以 及 塔里木 盆地 等 地 区 ,较为 显著 负 相 关 关 系 的 
区 域 ( -1.00 ~ -0.50) 主要 分 布 在 柯 坪 附 近 。 塔 
里 木 河 流域 NDVI 与 气温 呈正 相关 关系 的 区 域 占 全 
区 面积 的 27. 08% , 主要 集中 在 天 山中 部 .昆仑 山 
脉 .阿尔 金山 等 地 区 。 适 宜 的 气温 可 以 通过 光合 ZS 
腾 .呼吸 等 代谢 过 程 促进 植被 生长 变化 , 当 气 温 上 升 
超过 植被 的 适宜 温度 时 ,会 对 植被 生长 代谢 过 程 产 
生 抑 制作 用 ,从 而 NDVI 与 气温 呈现 负 相 关 。 塔 里 
木 河流 域 NDVI 与 气温 呈正 相关 关系 的 区 域 主要 分 
布 于 高 山区 ,这 可 能 由 于 高 海拔 地 区 气温 相对 较 低 ， 
热量 成 为 限制 植被 萌发 生长 的 主要 因素 。 
塔里木 河流 域 NDVI 在 大 部 分 区 域 (82.48% ) 
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7 2000—2018 年 塔里木 河流 域 NDVI 与 气温 、 
降水 量 相关 分 析 


Fig.7 Correlation analysis between NDVI and temperature 
and precipitation in the Tarim River Basin from 


2000 to 2018 


与 降水 量 呈 正 相关 关 系 ,几乎 覆盖 了 整个 区 域 , 与 降 
水 量 呈 正 相 关 关 系 较为 显著 的 区 域 (0.50 ~ 1. 00) 
主要 分 布 在 克 州 、 喀 什 地 区 。 塔 里 木 河流 域 NDVI 
与 降水 量 呈 负 相 关 关 系 的 区 域 主 要 集中 在 天 山中 
部 .昆仑 山脉 .阿尔 金山 等 山区 , 且 相 关 性 普遍 较 低 。 
水 在 植被 的 养分 传输 ,结构 维持 以 及 代谢 过 程 中 起 
着 十 分 重要 的 作用 ,同时 ,降水 增加 土壤 湿度 ,为 植 
被 提供 充足 水 分 ,从 而 促进 植被 生长 变化 。 对 高 山 
地 区 而 言 ,这 些 地 区 降水 量 较 高 ,对 植被 生长 相对 充 
足 ,温度 相对 较 低 ,植被 对 降水 量 的 响应 不 显著 。 
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图 8 2000—2018 年 塔里木 河流 域 NDVI 与 气温 、 


降水 量 偏 相关 分 析 


Fig.8 Partial correlation analysis between NDVI and 


temperature and precipitation in the Tarim River 


Basin from 2000 to 2018 


植被 的 变化 受到 诸多 气候 因子 的 综合 影响 , 单 
一 的 气候 因子 往往 不 能 解释 植被 复杂 的 变化 过 程 ， 
因此 引入 偏 相关 分 析 法 来 研究 气温 和 降水 量 对 植被 
的 共同 影响 。2000 一 2018 年 塔里木 河流 域 NDVI 和 
气温 的 偏 相 关系 数 ,与 NDVI 和 降水 量 的 偏 相 关系 
数 , 绝 大 多 数 (82.41% ) 在 -0.50 ~0.50 范围 内 , 相 
关 性 普遍 较 低 , 且 两 者 在 -0.50 ~0.50 之 间 的 空间 
分 布 差异 明显 。2000 一 2018 年 塔里木 河流 域 NDVI 
与 气温 和 降水 量 的 平均 偏 相关 系数 为 -0. 05 ,0.21。 
总 体 而 言 ,2000 一 2018 年 塔里木 河流 域 NDVI 与 降 
水 量 呈 现 正 相关 ,与 气温 呈现 负 相 关 , 晶 与 降水 量 关 
系 更 密切 。 进 一 步 说 明 , 对 干旱 半 干 旱 的 塔里木 河 
流域 而 言 , 水 分 在 一 定 程度 上 是 决定 植被 生长 变化 
的 主要 因素 。 
气温 .降水 量 与 植被 NDVI 的 偏 相关 性 的 高 低 
显示 了 气候 因子 对 植被 影响 力 的 强 弱 ,通过 对 偏 相 
关系 数 求 绝对 值 , 来 表示 气候 因子 对 植被 的 影响 力 
的 大 小 。2000 一 2018 年 塔里木 河流 域 气 温和 降水 
量 对 NDVI 的 综合 影响 力 从 东北 向 西南 逐渐 增加 
(图 9) ,气候 因子 对 NDVI 影响 力 在 0.00 ~0. 50 范 
围 内 居多 ,这 表明 气温 和 降水 量 对 NDVI 的 影响 相 
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9 2000—2018 年 塔里木 河流 域 气候 因子 
影响 力 分 布 ,NDVI 残 差 分 析 
Fig.9 The influence of climatic factors to vegetation 
and residual analysis of NDVI in the Tarim River 
Basin from 2000 to 2018 


对 较 弱 ,塔里木 河流 植被 的 生长 变化 还 应 与 其 他 相 
关 因 素 紧密 关联 ,如 塔里木 河流 域 复杂 地 质地 貌 .水 
系 分 布 以 及 人 类 活动 等 。 
2.4 塔里木 河流 域 人 类 活动 对 植被 变化 的 影响 

通过 气候 因子 对 植被 变化 的 综合 影响 进行 分 
离 ,得 到 2000—2018 年 塔里木 河流 域 NDVI 残 差 的 
空间 分 布 ( 图 9) 。 残 差 通 常 受 几 种 因素 影响 ,其 中 
地 形 等 固定 空间 因子 不 能 导致 残 差 趋势 的 定向 改 
变 , 只 有 长 期 的 人 类 活动 和 其 他 线性 变化 因素 才 可 
以 造成 残 差 的 趋势 性 变化 。 

2000—2018 年 塔里木 河流 域 NDVI JRX IE X. 
域 , 主 要 分 布 于 库尔勒 绿洲 、 渭 干 河 三 角 洲 、 阿 克 苏 
河 冲 击 平原 绿洲 .喀什 三 角 洲 绿洲 、 叶 尔 羌 河 绿洲 、 
和 田 绿洲 绿洲 范 漠 过 渡 带 以 及 塔里木 河 下 游 附 近 。 
NDVI 残 差 趋 于 正 ,表明 这 些 区 域 人 类 活动 等 其 他 线 
性 因素 对 植被 变化 产生 了 积极 的 影响 。 对 于 绿洲 
区 .绿洲 荒漠 交错 带 , 由 于 人 口 的 增长 ,人 类 不 断 开 
晨 土 地 ,使 得 人 工 绿洲 面积 不 断 扩大 ,从 而 促进 了 植 
被 覆盖 变化 。 对 塔里木 河 下 游 附近 而 言 , 自 2001 年 
起 ,塔里木 河 流域 生态 治理 工程 的 有 效 实施 ,在 一 定 
程度 上 保障 了 下 游 来 水 ,有 效 遏 制 了 塔克拉玛干 沙 
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漠 和 库 鲁 克 塔 格 沙漠 合拢 趋势 ,改善 了 7 下游 植 被 生 
境 , 从 而 使 得 植被 长 势 持 续 变 好 。 

2000—2018 年 塔里木 河流 域 NDVI 残 差 为 负 区 
域 , 主 要 分 布 于 哈密 地 区 .吐鲁番 地 区 EE 
沙漠 、 帕 米尔 高 原 、 昆 仑 山脉 .阿尔 金山 脉 等 地 区 , 残 
差 绝 对 值 相 对 较 小 ,这 些 区 域 人 类 活动 影响 相对 较 
弱 , 除 气温 降水 量 因 子 影响 外 , 则 可 能 由 于 其 他 持 
续 变 化 的 线性 因子 使 得 残 差 为 负 , 从 而 对 植被 生长 
变化 产生 消极 作用 。 对 于 人 口 稀少 的 亮 漠 地 区 而 
言 , 则 可 能 由 于 气温 持续 变 暖 加 剧 地 表 燕 散 过 程 ,使 
得 土壤 干 层 发 育 ,对 植被 生长 变化 产生 抑制 作用 。 
对 于 高 山区 等 山体 多 地 带 , 则 可 能 由 于 其 气候 条 件 、 
地 质 环 境 等 因素 (山区 常年 积 雪 、 恶 劣 地 质 条 件 
等 ) ,对 植被 生存 、 生 长 具有 消极 的 影响 。 


3 结 论 


本 文 基于 MODIS-NDVI 数据 ,利用 Sen + Mann- 
Kendall 趋势 分 析 、 变 异 系数 法 、Hurst 指数 、 偏 相关 
分 析 及 残 差 分 析 等 方法 ,对 2000—2018 年 塔里木 河 
流域 植被 覆盖 时 空 变化 特征 、 植 被 可 持续 性 进行 分 
析 , 结 合同 期 气候 数据 ,探究 气候 因子 和 人 类 活动 等 
因素 对 塔里木 河流 域 植被 变化 的 影响 ,得 到 以 下 
结论 : 

(1) 从 时 间 上 看 ,2000 一 2018 年 塔里木 河流 域 
植被 覆盖 总 体 呈 现 明显 的 增加 趋势 。 植 被 NDVI 以 
2008 年 为 转折 点 ,存在 明显 的 阶段 变化 。 从 空间 上 
看 ,塔里木 河流 域 植被 分 布 具 有 明显 的 地 域 差 异 。 
“ 北 高 南 低 , 西 高 东 低 ”, 植 被 覆盖 主要 分 布 于 天 山 
山脉 、 帕 米尔 高 原 以 及 昆仑 山脉 东部 绿洲 地 区 及 绿 


洲 荒 漠 交 错 带 。 男 外 ,流域 内 植被 沿 河 流 沿岸 分 布 ， 
呈现 出 条 带 状 。 


(2) 从 NDVI 变化 趋势 上 看 ,塔里木 河流 域 大 
部 分 区 域 植被 覆盖 趋势 基本 不 变 , 植 被 改善 的 区 域 
占 全 区 面积 的 37. 97% ,退化 区 域 占 11.24% ,山区 
水 资源 丰富 ,加 之 气温 变 暖 ,山区 植被 活动 响应 
显著 。 

G) 从 植被 覆盖 可 持续 性 分 析 上 看 ,塔里木 河 
流域 大 部 分 植被 覆盖 趋势 基本 不 变 。 反 持续 区 域 比 
重 为 30.83% , 主要 分 布 在 山区 地 带 和 水 系 附近 ; 持 
续 性 区 域 比重 为 69.17% ,多 分 布 在 绿洲 .绿洲 荒漠 
交错 带 以 及 荒漠 区 。 

(4) 塔里木 河流 域 气温 和 降水 量 对 NDVI 的 综 


合影 响 从 东北 向 西南 逐渐 增强 。 总 体 上 ,NDVI 与 降 
水 量 呈 正 相 关 性 ,与 气温 呈 负 相关 性 , 且 与 降水 量 的 
相关 性 更 强 , 降 水 量 对 NDVI 的 影响 更 为 明显 。 人 
类 活动 对 植被 变化 产生 积极 影响 的 区 域 主要 分 布 于 
绿洲 ,绿洲 范 漠 过 渡 带 .塔里木 河 下游 附 近 。 芹 漠 以 
及 山区 地 带 等 人 类 活动 相对 较 弱 的 地 区 ,其 NDVI 
残 差 还 受 其 他 持续 变化 的 线性 因素 影响 。 本 文 仅 探 
讨 气温 .降水 量 气候 因子 以 及 人 类 活动 对 塔里木 河 
流域 NDVI 的 影响 ,由 其 残 差 分 析 结 果 可 知 ,塔里木 
河流 域 NDVI 的 变化 还 受 其 他 空间 线性 变化 因素 影 
响 ,如 土壤 湿度 .水系 径流 、 积 雪 获 羡 等 其 他 因素 ,后 
续 将 进一步 研究 塔里木 河流 域 NDVI 的 其 他 影响 
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Spatial-temporal variations of vegetation coverage in 
the Tarim River Basin from 2000 to 2018 
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Abstract: The ecology of the Tarim River Basin, Xinjiang, China is relatively fragile. During the past few dec- 
ades ,unreasonable development and utilization have resulted in serious damage to the vegetation coverage. Due to 
the implementation of the ecological comprehensive management project in the Tarim River Basin and by recent cli- 
mate change ,the temporal and spatial distribution of vegetation in the basin has undergone certain changes. Clarif- 
ying the relationship between vegetation coverage change and climate change of the basin, as well as human activi- 
ties,can provide a scientific reference for ecological maintenance and management in the Tarim River Basin. 
Hence , based on MODIS-NDVI and using the Sen + Mann-Kendall trend analysis, coefficient of variation method, 
Hurst index, partial correlation analysis, and residual analysis, this study analyzed the temporal and spatial change 
characteristics and sustainability of vegetation coverage in the Tarim River Basin from 2000 to 2018 and investigated 
and distinguished the range and extent of the impacts of temperature and precipitation climatic factors as well as hu- 
man activities on vegetation coverage changes. The results are summarized as follows: (1) The level of vegetation 
coverage in the Tarim River Basin has increased significantly with clear stage changes. The vegetation changes var- 
ied from an increasing trend to a decline before 2008 and then increased rapidly. The vegetation increase is mainly 
concentrated in Korla, Yangi, Aksu, Kuga, Kashi, and the region of the central Tianshan Mountains. The vegetation 
variation near the Kuluktag Mountains is relatively small. (2) More strongly affected by regional geomorphologic 
patterns ,the vegetation coverage of the Tarim River Basin in the north is greater than that in the south, and the veg- 
etation coverage in the west is greater than in the east. The area of relatively high vegetation coverage is distributed 
in the mountainous areas, oasis areas, and the oasis-desert interlaced zone. In addition, the basin vegetation is dis- 
tributed along the river, shaped as a strip. (3 ) The vegetation coverage trend in the Tarim River Basin is mostly un- 
changed (50.80% ) and the improved trend area (37.97% ) is mainly concentrated in the oasis with abundant wa- 
ter resources in the basin. The stability of vegetation coverage in the Tarim Basin is in a state that “the outer side of 
the basin is higher than the inner side and high and low fluctuations coexist and the moderate fluctuation is most- 
ly”. The response of vegetation activity is significant in areas covered with high vegetation amounts such as moun- 
tainous areas. Most regions in the study area will maintain the existing stable trend. The area of unsustainable vege- 
tation coverage change (30.83% ) is mainly distributed in mountainous areas and water systems and the sustainable 
area (69.17% ) is mostly concentrated in oasis , oasis-desert interlaced zone ,and desert area. (4) The NDVI in the 
Tarim River Basin has a significant negative correlation with temperature in most areas (72.92% ). The NDVI is 
positively correlated with precipitation in most areas (82.48% ).The combined effects of temperature and precipita- 
tion on the NDVI gradually increase from northeast to southwest and the impact of precipitation on the NDVI is more 
significant. The areas where human activities have positive effects on vegetation change are mainly distributed in the 
oasis , oasis-desert transition zone, and the lower reaches of the Tarim River. In areas with relatively low human ac- 
tivity such as deserts and mountainous areas ,the NDVI residuals are also affected by other linear factors of continu- 
ous change. 

Key words: vegetation coverage; spatial-temporal variations; climatic factor; human activity; the Tarim River 
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